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ABSTRACT 


A new method for absolute calibration of electroacoustic trans- 
ducers is described which permits low frequency reciprocity calibration 
in a test apparatus of small dimensions. The method takes advantage of 
the reduced wave speed in a cylindrical water column bounded by a 
compliant PVC tube to reduce the length of the standing wave to 
dimensions which can be easily handled in the laboratory. 

Results of a limited series of experiments show a mean reproduci- 
bility of +0.8 decibels in a one-meter long tube at frequencies between 
750 and 1100 Hz and sound speeds from 330 to 365 meters/second. The 
hydrophone voltage sensitivities obtained by the method, however, 
systematically differ by 1.5 decibels from those measured by comparison 
with a standard hydrophone; additional research will be required to 


resolve this discrepancy. 
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1. WHY COMPLIANT TUBE RECIPROCITY WORKS 


A. ACOUSTIC CALIBRATION AND RECIPROCITY 


There are two general classes of underwater acoustic transducer 
calibration: absolute, or primary, calibration, and comparison, or 
secondary, calibration. The latter presumes the availability of a sound 
field of Known amplitude measured by a calibrated transducer, but 
absolute calibration can be performed using only voltage measurements 
and primary measurements of the test system environment in a sound 
field of unknown amplitude. Once the voltage sensitivity of a standard 
hydrophone has been established in a particular frequency range and 
environment by absolute calibration, it may be used as a reference for 
comparison calibration of other hydrophones and of projectors; this 


absolute calibration is most often done by the reciprocity method. 
1. Fundamentals of Reciprocity Calibration 


The reciprocity method for absolute calibration is based on the 
properties of linear, passive systems whose action can be described by 
the two-port equation (1.1), which relates the KS Tie pressure and 
velocity (p,u) at a transducer's mechanical terminal and the voltage and 
current (E,l) at its electrical terminals by a symmetric or antisymmetric 


matrix of constant coefficients; such a system is called reciprocal. 


PLIN . 12 m 


Hunt [Ref. 1] derives this equation from first principles and shows how 
the matrix elements may be evaluated for common transducers. Foldy 
and Primakoff show [Ref. 2] that the short circuit transmitting current 
response S and the open circuit hydrophone voltage response M of a 


reciprocal transducer T are related by a constant factor 


(1.2 8 


s 
ዘ 
SUE 
ዘ 

"U |= 


T 


where J is called the reciprocity parameter and is characteristic of the 
test medium and its boundaries. Bobber has demonstrated [Ref. 3] 
that J is eguivalent to the acoustic transfer admittance between an 
arbitrarily defined transmitting surface T generating a volume velocity 
Ú and an arbitrarily defined receiving surface H within the test envi- 
ronment exposed to an average pressure P, and exists as a function of 
the test medium and its boundaries for any arbitrary linear, passive, 
and reversible medium. 

| When two test transducers? are located at points im WNE 
medium where the pressure amplitudes are equal, the reciprocity param- 
eter allows expression of the signal voltage of a test hydrophone H 
placed at one point as a function of the transmitter current response of 
a reversible transducer T placed at the other point and driven to 


create the sound field; 


Етн = МыРы = MuSrlr (1999 


where | is the drive current input to T. If the sound field generated 
by T is then replaced by a field which is identical at the locations of T 
and H, but is generated by an separate projector P, the outputs of T 
and H in this new field will be related by equation (1.4). Combining 
this equation with (1.3) produces the calibration equation (1.5), which 
allows the computation of the test hydrophone's voltage sensitivity H 
directly from voltage measurements made in the three source-to-receiver 


configurations:  T-to-H, P-to-H, and P-to-T. 


*To avoid confusion, the term test transducers will always be used 
when referring to both transducers used in the reciprocity, calibration. 
The term reversible transducer, or the ANSI designator T , will be 
used to refer to the reciprocal transducer used both as a hydrophone 
and as a projector during the calibration; test hydrophone, or H will 
designate the transducer used only as a hydrophone; and projector, or 
P will identify a third sound source used to create the sound field 
used in the test. 
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Е. МР Ј5, (1.4) 
Een MP -My 
| (1.5) 
MA = | E pu Ery 
Eor ly 


2. Reciprocity Test Geometries 


Most reciprocity calibrations are made in one of the two 
earliest system geometries for which the reciprocity parameter was 
computed: free field [Ref. 4], ar close-coupled [Ref. 5]. The tech- 
niques for carrying out these tests, and the limitations required to 
ensure the validity of the test are described by the ANSI standard for 
underwater calibration [Ref. 6]. The free field test geometry requires 
a volume of water large enough to permit spherical spreading of the 
acoustic waves undisturbed by boundary reflections. This implies tanks 
of dimensions greater than ten wavelengths (about 15 meters at 1 kilo- 
hertz) for underwater testing. Coupler reciprocity is limited by the 
requirement to encase the entire active surface of two test transducers 
in a rigid-walled cell of dimensions less than 1/10 of a wavelength; 
useful underwater low frequency transducers limited to such dimensions 
2507 Ce difficult to build. 

The practice of reciprocity calibration is not, however, 
constrained to the two conventional geometries. Rudnick [Ref. 7] 
discusses three compact geometries with known reciprocity factors; 
Bobber's reference text on measurements [Ref. 8] describes several 
more. These unusual test environments extend the practice of reci- 
procity calibration beyond the limits of of the two conventional geome- 
tries. A summary of these methods is included in Table |; each of the 
reciprocity parameters listed is the quotient of an area characteristic of 
the system geometry, and of the specific acoustic impedance of the 
medium. 

The compliant tube test system described in this paper adds 
another option to the array already available to the experimenter. It 


has particular value to the underwater acoustician because it 


EE 


Some Reciprocity Calibration Systems 


Geometry Reciprocity Parameter Reference 
Free field ” 2dA 4 
Poc 
Close—coupled € 2 
. Hat 
Plane-wave — 8 
PoC 
Cylindrical -wave — 8 
PoC 
Helmholtz resonator 2nALi 7 
PoC Q 
Waveguide, TALT 7 
resonant plane wave poc Q 


significantly lowers the frequency limit for underwater reciprocity cali- 
bration below the practical limits of the conventional geometries. 
Unlike rigid-tube methods, it can be accomplished within a convenient, 
Inexpensive test cell, and although still limited by the requirement for 
accurate measurement of the system quality factor, does escape resonant 
reciprocity calibration s usual restriction to discrete frequencies. The 
“virtual rigid walls” created by the standing wave pressure antinodes, 
and in which the transducers are placed for the test, can be adjusted 
to any resonant length and frequency by simply raising or lowering the 
water level. 

Our development of the compliant tube reciprocity parameter 


must begin from an understanding of wave motion in the water column, 


and of the oscillatory expansion of the tubes circumference. From 
these equations, we will develop the energy balance of the oscillating 
system; this in turn will provide the result for the reciprocity 


parameter. 


B. COMPLIANT TUBE WAVE EOUATION 
1. Waves in the Water 


Lighthill's development [Ref. 9] of the wave equation for a 
lossless fluid moving coaxially in a lossless tube of varying cross 
section includes the case of the compliant elastic tube in which our 
experiment is performed. The momentum and mass conservation equa- 


tions for an inviscid, ideal fluid are 


2 2) = -2 (1.6) 
öt öz öz 
(pa) WIE 7) 





T + — (pau) = 0 
where p is the acoustic excess pressure, u the acoustic coaxial particle 
velocity, and a the local time-varying cross-sectional area of the tube, 
which varies slightly about its equilibrium value A. When fluid entropy 
is constant, both the local density p and the local cross section a may 
be expressed as functions of the acoustic pressure. (The latter relation 
is just the inverse of the function which relates pressure to cross 
section variations.) Discarding terms of second order and higher 


reduces these eguations to 


õu Әр (158) 
55 

(pa) _ „Bu (1.9) 
T 7: 
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Taking the spatial derivative of (1.8) and the time derivative of (1.9), 
combining the results with the linearized equation of state p = (p -p,)c* 


and discarding terms above first order gives a wave equation 


ፅፀ -2 0%p 
መመመ == с . 
8z a" (1.105 


where the wave speed c is defined by 
: 13 
K pol CS S Í A p = ol k A d (1:8 
Geh Op /, Poópl, A õp], 


(The subscript "0" above implies evaluation of the derivative at the 
undeformed wall position and under adiabatic conditions.) The first 
term of equation (1.11) is the adiabatic compressibility or reciprocal 
bulk modulus K; if there were no dependence on cross section, this 
would reduce to the equation for the sound speed in the bulk medium. 
If the second term, called the distensibility D, is much larger than the 
adiabatic compressibility, the speed of sound waves in the tube will be 
determined by the characteristics of the wall rather than by those of 
the bulk fluid, and will be much slower. This in turn will reduce the 
size of the resonant tube reguired to operate at a given freguency. 

The harmonic solution of eguation (1.10) with a pressure 
release boundary at z=0 and a rigid boundary at z=L, and for frequen- 


cies below the cutoffs for azimuthal and radial modes, is 
p(z,t) = P sin(kz)sin(wt) (1.12) 


with a corresponding velocity field 
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p 
u(z,t) Ze cos(kz) cos(wt) (1 13) 


0 


The wave number k is limited by the boundary conditions to the 


discrete values 


em eut MI 023... | (1.14) 


ke = 
A L 


This sound field has associated with it both kinetic and potential energy 
density. The kinetic energy of axial water motion in a thin radial cross 


section of length dz is 
pa” 
da = 57 (1.15) 


The potential energy in the same section is found by considering the 


work done by the fluid as the tube boundaries expand: 


Ç ና 1 3 
dE, = — | pav = -| ንይ መባ -adz E ሠ GNIS 
0 0 ል 951. .. ÖP, 


2 GE 
= (K - 0) Adz E 


where ž is a coordinate of the periodic motion. 


2. Motion of the Wall 





The distensibility can be calculated from the material proper- 
ties of the tube wall and its geometry. In the case of a cylindrical tube 
which responds elastically to small excess stresses, the circumferential 
tensioning force component due to the pressure on a thin radial cross 


section of length dz is 


is 


1/2 
F = 2Rdz | p cos 8 d8 = 4Rdz p 
—n/2 


(1-15) 


00.55. ሎ. (1.18 


6 A 2HYdz BY 


This tensioning force causes a strain e defined by equation (1.19). 
Where p is small and a varies only slightly from.its undisturbed value 
A, the distensibility is then 


D 


! da 
d 


A dp hY 


| 55 (1.20) 

0 

This stretching of the vvall along its circumference has the characteris- 
tics of a damped, stiffness-controlled driven oscillator (neglecting long- 
itudinal effects). The ratio of the driving pressure at any cross 
section and the circumferential strain velocity is the mechanical impe- 
dance of the cross-sectional ring treated as a damped oscillator of 


negligible mass with elastic constant Y and damping coefficient R(m) 
(Ref. 10:pp.10-12] 


ne 2P sin(kz)sin(wt—@) 


A | (1.21) 
~N 2 2 - (1 . 220 
W = x İRA “(Y/cə) ያሙ ureisi əl = 
cd m 


The polyvinyl chloride (PVC) tube used in this experiment is actually 
not elastic over any finite range of distortion, so the parameter Y is 
not Young's modulus; it is only a linear approximation to the PIaSNIE 
deformation that the PVC undergoes in harmonic motion. Within this 
linearization, the potential energy stored in a thin radial section of the 
wall as the tube is expanded by a small amount from its initial cross 


section can be found from the force-displacement integral 


t 
dEw - | ፍወ d a ор” Adz (22) 
0 Ç 


with the force given by equation (1.18). 


7: COMPLIANT TUBE RECIPROCITY PARAMETER 


The arguments used by Garrett, et al., [Ref. 11] to derive the 
reciprocity parameter J for a rigid tube will also yield an expression for 
the reciprocity parameter of the free surface, compliant wall tube. As 
with the rigid tube, the solution depends on relating the energy 
content of the system E and the energy loss per cycle AE to the 
acoustic transfer admittance (reciprocity parameter) described by equa- 
tion (1.2). When the ratio of energy content to energy loss rate is 
known in terms of 2 and of primary measurements, the resonance 
quality factor equation (1.24) (which can also be found from primary 


measurements) can be used to evaluate J. 








ን ድሜ (1.24) 
AE 27 


1. System Energy Content and Losses 


The compliant tube system stores energy both in the water 
column and in the tube wall. Kinetic energy of motion of the water 
pressurizes the water locally, which in turn stresses and expands the 
tube diameter. The energy balance exchanges between the three 
degrees of freedom are governed by the work equations for pressuriza- 
tion (equation (1.17)) and for expansion (equation (1.23)), and by the 
kinetic energy equation (1.15). Summing these equations gives the 


total energy in a thin cross section 


2 2 2 
dE = dE, + dE, * dE, = A dz | Pall (k - 0)5+ Dp (1:25) 


2 2 u” ER (1.26) 
z pat LE + o) b= adz ) Pa + „af 
o dz 2 T ) 2 209€ 


where equation (1.11) has been used to replace K and D with p,c?. 
The sound field in the absence of attenuation is described by equations 
(1.12) and (1.13); substitution gives a form which can be integrated 


over the length of the water column 


L 2 
E = | iə () cos (kz)cos (ut) E 2 sin (kz) sin^(ot)? dz. (1.27) 
ር 2095€ | 


ALP * (1.28) 
ə 


This is identical to the result for a rigid resonant tube [Ref. 11]. 

In the steady state, the cyclic energy loss equals the work 
done by the driving transducer. If this transducer is located at a 
pressure antinode, pressure and velocity must be in phase at the 
transducer face ( since the impedance to flow at a velocity node is infi- 
nite). The energy. transferred to the system is the work done, which 


can be written as 


ə Do us A ds =P Usin (vt) (1.29) 


where U is the volume velocity amplitude of the transducer surface. 


Integrating over one cycle gives 


PU (1.30) 
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2. Reciprocity Parameter 


Our calculations above have yielded both the system energy 
content and the system energy loss (actually, its throughput at steady 
state). The energy equation for Q (1.24) immediately produces 


E "PAM 


AE U 2pyc” 





= 61.31) 
27 

The ratio of volume velocity amplitude to pressure amplitude that 
appears in this result is the acoustic transfer admittance between pres- 
sure maxima (equation (1.2)), and must also be the reciprocity param- 
eter for hydrophones positioned at separate maxima when the column is 
driven at resonance. Equation (1.32) follows; again, it is the same as 


the reciprocity parameter for the rigid tube. 


ję TALf 
poe a ECH 


The solution of equation (1.32) requires the experimental measurement 
of three resonance parameters: the resonance frequency f; the reso- 
nance quality factor Q; and the wavelength' of the resonance standing 
wave, in order to compute the sound speed C. The accuracy of the 
value computed for J depends on our ability to measure these parame- 
ters accurately. The methods for doing this will be discussed in the 
next chapter. 

In using eguation (1.5) with the result derived above, we 
must ensure that the pressures and voltages implied are indeed the 
maximum amplitude values at the antinodes. The test transducers: 
output will be proportional to the average sound field over their active 
surfaces, which at its maximum will be less than that which would be 
generated by a uniform field at the maximum pressure amplitude P by 
the factor 
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_ sin(kL, /2) 


kL./2 (1.33) 


1 
71 =— $ cos k(z-zm) ds 
Ar 


Measured voltage outputs must be corrected by this factor before the 
results are used in the calibration equation when transducer dimensions 
are an appreciable fraction of a wavelength. For the LC-10's used in 
our experiments, this factor is 0.05 dB at 1kHz; for the LC-32's, it is 
02 siis 


3. Limits of Derivation 


The derivation of the compliant tube reciprocity parameter has 
been done under four important simplifying assumptions. To the degree 
that these approximations do not represent the real system, the results 
obtained by absolute calibration will be in error. 

The first of these approximations is the use of the solutions 
to the lossless wave equations in computing the total energy in equation 
(1.28). The PVC wall does attenuate the sound wave, resulting in a 
resonance sound field which resembles that derived by Kinsler, Frey, 
Coppens, and Sanders [Ref. 10:pp.206-209]. Figure 1.1 shows the 
measured variation of the pressure amplitude with depth for the n=4 
mode. The pressure amplitudes at separate maxima will therefore not 
be identical, as is required by our development, which assumes an infi- 
nite standing wave ratio. 

A second approximation is the assumption of linear behavior 
for the polyvinyl chloride wall. In fact, PVC has no elastic region for 
strain, as was discussed in the derivation of the wall motion equation. 
The nonlinearity of the wall permits mode coupling into other resonance 
modes, particularly when two resonances have spatial pressure maxima 
near the same position and are related subharmonically. The result of 
such coupling when the column is driven at one frequency is to depress 
the amplitude of individual pressure maxima by an amount proportional 
to the degree of coupling at that position; again violating the require- 


ment for equal pressure amplitudes at the two test positions. 


Third, the reciprocity concept requires that the receiving 
transducer approximate an infinite mechanical impedance to the medium. 
While piezoceramics do meet this approximation very well, the soft 
mounting material used the the hydrophone body and gas bubbles on 
the hydrophone surface do not. Their presence could significantly 
distort the sound field in the vicinity of the transducer. 

Finally, the projector used to drive the column into resonance 
does not produce a sound field identical to that generated by a small 
piezoelectric transducer suspended in the water column. The local 
sound field near the cylinder is some combination of radial modes which 
match the cylindrical boundary conditions to the overall plane wave, 
coupling the driving transducer to the standing wave system. The 
effects of this coupling factor are slightly frequency-dependent, and 
are not accounted for in this model. 

Our measurements demonstrate that the simple model derived, 
despite these deficiencies, produces consistent results to within +0.8 
dB, but does not explain differences from a standard comparison cali- 
bration of up to 2.1 dB. The root-mean-square average difference is 
1.5:0.4 dB for the range of hydrophone types and frequencies tested. 
A method and apparatus for performing this type of calibration experi- 


ment, and the results of our experiments, follow. 
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Figure 1.1 Pressure Amplitude Depth Profile (n=4) 
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11. COMPLIANT TUBE NECIPROGITYFENFERIMENTS 


A. RESONANCE PARAMETER MEASUREMENT 


The formula for the reciprocity parameter derived in the preceding 
section requires the accurate measurement of three resonance parame- 
ters: the resonance frequency, the sound speed at that frequency, 
and the resonance quality factor: Each of these measurements involves 
finding an extremum from a set of .data samples, which will inevitably 
contain noise. The experimenter who attempts to resolve extrema from a 
discrete sample set by comparing measured signal outputs performs the 
discrete equivalent of differentiation, which enhances noise and 
increases net error in the result. 

The solution to accurate, reproducible resonance parameter meas- 
urement which D. Conte and S. Garrett originally developed [Ref. 12] 
combines the high data rate of automated measurements with the power 
of curve fitting. The intuitive explanation of the effectiveness of this 
method is the human analogue; rather than take discrete measurements, 
experimenters instead "swing the system slowly through the 
extremum, and let the eye interpolate the location of that point. The 
optimum automated processing analogue for this method is curve fitting 
to the minimum expansion of any function y(x) which will resolve 


extrema, the second order equation 
y(x) = Co dE Cix T CZE (2 1) 


The three coefficients of this curve can be found for any set of three 
or more measurements; if more than three measurements are made, the 
nonlinear regression methods described by Bevington [Ref. 13] can be 
used to reduce the noise error in the measurement of the extremum, as 


long as the span of the sample set around the extremum is small enough 


that equation (2.1) is a valid approximation. The regression equations 
for the three coefficients are algebraic, and are derived in Appendix B 
(equations (B.12), (B.13), and (B.14)). The location of the extremum 
(again from Appendix B) is displaced from the coordinate origin of the 


sample set by 
Xu 00 -C,/2Cz (2.2 


İn our application of this method to measurement of the resonance peak 
frequency , the sample set is composed of seven hydrophone output 
voltage amplitudes measured at equal frequency intervals; the sample 
set origin is shifted until the estimated peak frequency falls within one 
interval of the center frequency of the sample set. This is the peak 
frequency locating algorithm. Figure 2.1 is a example of this technique; 
the asterisks mark data points used to establish the coefficients for the 
curve and to compute the location of the maximum, and points marked 
with circles are discarded from the final fit because they are more than 
three steps from the extremum. The sound speed is found by 
sampling the pressure amplitude at regular intervals in depth from the 
surface, and again shifting the set used for curve fit until the esti- 
mated maximum or minimum is within one step of the center of the set. 
The positions thus located occur at quarter wavelength steps, and can 
be fit to 


2 = MAN/4 + Zo ሠ... 


to obtain the wavelength; the sound speed is then just the product of 
the resonance peak freguency and the wavelength. This is the sound 
speed measurement algorithm. 

The curve fit in freguency also yields a reduced-noise value for 


the peak amplitude at resonance, which is essential to accurate 


measurement of the half-power frequencies required for the Q measure- 
ment algorithm (equation (1.24)). Once the half-power amplitude is 
known, the half-power frequencies can be found by sampling outward 
from the peak frequency until a pair of samples brackets this value at 
the minimum frequency interval permitted by the equipment (in our 
case, 1 Hz). The half power frequency is then calculated by linear 
interpolation. 

The curve-fitting method can resolve extrema in noise to a degree 
far surpassing the eye integration method it was drawn from. İt has 
the additional advantages of sound mathematical basis, speed, reprodu- 
cibility, and independence from observer biases. It also has the consid- 
erable advantage of being easier on the observer. The method's main 
drawback is its bias toward finding the parabolic curve in any data, 
including pure noise; care must be taken in selecting stepsizes to 
ensure that the changes in the signal are not so small that they are 
overwhelmed by noise, and to ensure that the sample set span is not 
larger than the region about the extremum where a second order 
approximation is valid. The best method for checking new cases is 
graphical display of the data set and the curve for the operator 


approval. 


DESCRIPTION OF EQUIPMENT 


1. Control System (Figure 2.3) 


The experiment is controlled by a venerable Hewlett Packard 
HP-85 laboratory computer with 32K memory extension, !EEE-488 bus 
controller, and appropriate ROM extensions. The program used is 
documented in Appendix C. Commands generated by the program over 
the data bus control the experiments configuration changes during 
measurements through the two dual remote coaxial switches (Hewlett 
Packard HP 59307) and a Hewlett Packard HP 3421 Data Acquisition 
Control Unit. The latter unit is configured with two actuator relays 
and eight measurement channels (option 020). The two actuator relays 


control the operation of the reversing DC manipulator drive motor. 
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Figure 2.3 
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Transmitter/ Receiver Schematic 


Figure 2.4 


The manipulator's screw drive, which is described below, also 
drives a Hewlett Packard HED 680 position encoder that increments a 
Hewlett Packard HP 5016 universal counter in totalizing mode. The 
changes in this totalized count provide (via the data bus) the program 


with a measurement of actual manipulator motion. 
2. Transmitter (Figure 2.4) 


A Hewlett Packard HP 3314 function generator provides a 
stable bus-controlled sine input signal at two voltage/current levels, 
and a synchronization signal at the selected frequency and constant 
voltage. The high voltage output is amplified and used to drive the 
piezoceramic transducer at rms voltages of 10 to 20 volts AC. The 
low-voltage, high-current demands of the USRD J11 projector [Ref. 14] 
are met by low-gain power amplification of the function generators low 
voltage output. The transmit remote coaxial switch (714)? selects the 
signal source from these two options and feeds the signal to the 
reversible transducer switch (706A), to the J11 projector, or through a 
test attenuator which returns a known signal to the receiver input 
switch (7068). 

A Hewlett Packard HP3478 digital multimeter, configured as a 
remotely controlled 5 1/2 digit AC voltmeter, is connected to the 


transmit signal line to measure the transmitting voltage. 
3. Receiver (Figure 2.4) 


The test transducers used in these experiments are: two 
Celesco LC-10 piezoceramic hydrophones (Figure 2.5), selected for their 
small size and good transmitting response; and a Celesco LC-32 piezo- 
ceramic hydrophone (Figure 2.6) of larger size. The voltage sensitivi- 
ties of all three units were measured by comparison with a reference 
standard USRD H56 hydrophone [Ref. 14]. The results are listed in 
Table 11. The LC-10's are mounted on 0.375 inch thinwall copper 
tubing; the LC-32 is mounted on 0.4375 inch tubing. Signal cables are 


?The remote switches are identified here by their IEEE-488 bus 
address codes. 
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TABLE Il 


Comparison Calibration Resutts 


(All values + 0.8 decibels) 


Frequency LC-1022166 LC-10#2338 1С-32285 
(Hz) 6 2m cable & 9m cable & Ilm cable 
740 -206.5 -206.8 2019-3 
860 -206.8 -207.3 -213.4 
1000 -206.9 -208.8 -214.3 


routed up the interior of the tubing to minimize effects on the sound 
field. The two test transducers are suspended from the manipulator 
platform. Their output is returned from the top of the manipulator 
assembly by two of three coaxial cables of the same length and meas- 
ured capacitance. The third of these cables connects the test attenuator 
into the receiver input switch (706B), which makes all receive signal 
path wiring nearly equivalent. The reversible transducer switch (/06A) 
selects either its receive or transmit lines; the receiver input switch 
(706B) selects the lockin analyzer input from the reversible transducer 
receive line, the test hydrophone, or the test attenuator. 

The EG&G PAR model 5204 lockin analyzer filters out the 
components of the received signal which are in-phase and in-quadrature 
with the reference (synchronization) signal. Their Pythagorean sum is 
used to obtain the magnitude. The analyzer is typically operated with 
an integration time of 300 milliseconds, which corresponds to an equiva- 
lent noise bandwidth of 0.8 Hz; this causes the lockin analyzer to act 
as a very sharp filter synchronized to the transmitter signal. The 
analog output of the model 5204 is passed to a measurement channel of 
the HP 3421 for digitization. 


4. Mechanical System 


The tests are conducted in a 6-inch (154 mm) inside diameter 
Schedule 40 polyvinyl chloride tube mated to the active surface of a 
USRD J11 electrodynamic underwater projector [Ref. 14]. The tube 
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. Figure 2,5  Celescóo EC-IUPRZATOPNENE 
(Reprinted with permission of Celesco Industries, Inc.) 


stands 1180 mm above the J11's active surface, and is filled to a depth 
of about 1 meter; its wall, which is smooth on the interior except for 
small regular extrusion made is 8 mm thick. The configura n 
resembles the USRD G40 calibrator used for comparison calibrations 
[Ref. 14]. (This system was originally constructed by Mills and 
Garrett for their work [Ref. 15] with fiber optic gradient hydro- 
phones.) Equalization for the J11 driver is provided by air pressuriza- 
tion of the installed air bag, monitored with a water-filled U-tube 
manometer. 

The test transducers are suspended on thinwall tubing 
supported on a rigid table which is raised and lowered on two one-inch 
rod tracks by a motor-driven screw drive; maximum range of motion is 
800 mm. The drive controls were described in the paragraphs on the 


control system above. 
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= Figure 2.6  Celesco LC-32 Hydrophone 
(Reprinted with permission of Celesco Industries, Inc.) 


MN? ROCEDURES FOR EXPERIMENT 
1. Preparation 


Before the test transducers are suspended in the water 
column, the experimenter must measure the reversible transducer s 
input capacitance, so that the transducer current can be computed from 
the transmitter input voltage and drive freguency. The active surface 
lengths of both hydrophones should also be measured for use in 
computing the transducer amplitude averaging factor (equation (1.33)). 
The test transducers are then suspended from the manipulator, wetted 
with a detergent solution, immersed in the water column, and brushed 
to remove bubbles from the active surface. If necessary, the tube's 
inner wall should also be brushed free of bubbles. The water column 


depth is measured with a meter stick. 


oo 


2. Mode Selection 


Two parameters are available to vary the resonant frequency 
of the water column at which the calibration is done: the water depth 
L and the mode number n (equation (2.4)). The procedure requires 
wavelengths much larger than the coaxial lengths of the transducers. 
In particular, the wavelength should be much longer than the coaxial 
length of the reversible transducer T to ensure that T's surface is 
entirely within the velocity node and feels little or no impressed pres- 
sure when T is transmitting. The vvavelength should also be at least 
four times the length of the test hydrophone H. Practical experience 
has shown that calibration results are most accurate when T and H are 
positioned at two different pressure maxima of the n=3 or higher order 


resonance modes. 


L = + 1/2) == (n + 1/2) (2.4) 


> Resonance Parameter Measurements 


After positioning the test transducers at the approximate loca- 
tions of the deepest and next deepest accessible spatial pressure maxima 
with similar pressure amplitudes, the operator uses the peak frequency 
measurement algorithm to find the resonance frequency f. The trans- 
ducers are then withdrawn above the first pressure minimum. 

One of the transducers is now used to measure the locations 
of maxima and minima in the standing wave pattern. Linear regression 
of this data gives the guarter wavelength, which in turn gives the 
sound speed. 

Finally, the transducers are precisely repositioned at the 
maxima where the resonance freguency measurement was made. A 
precise measurement of the peak frequency and of the Q factor are 
made, completing the set of measurements needed to compute the reci- 


procity parameter from equation (1.32). 
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4. Calibration Measurements 


Before making the final measurements, the experimenter 
switches the attenuated test signal into the receiver output to measure 
the net effects of cable losses, preamplifier gain drift, and deviations 
between the lockin analyzer's and the digital voltmeter's measurements. 
This correction is applied along with the amplitude averaging correction 
factors from equation (1.33) to the calibration voltage measurements 
that follow. 

The transmitter is switched to the reversible transducer T, 
and the receiver to the test hydrophone H. The transducer-to- 
hydrophone output voltage E(TH) is measured, corrected, and 
recorded, as is the transmit voltage V(T). The latter is used vvith the 
capacitive admittance of T to compute the transmitter current 1(T). 

The transmitter is then switched to the projector, and the 
drive voltage adjusted to give a sound field of the same amplitude at H 
as was used in the previous step. The projector-to-hydrophone output 
voltage is recorded (it is equal to E(TH), within the limits of the 
projector drive adjustment error). The final calibration measurement is 
made with the receiver switched to T with the projector drive 
unchanged; the result is the projector-to-transducer output voltage 
E(PT). | 

These measurements complete the reguirements to compute the 
test hydrophone voltage sensitivity from eguation (1.5). The same set 
of measurements can be used to compute the reversible transducer's 


sensitivity using 


2 E E 


— PTIB 
R (2.5) 


al 


The results for M are in volts per pascal; the customary units quoted 


are decibels re 1V/uPa, which is found by calculating 


Al = 20 log M — 120 (2.6) 


39 


5. Self-consistency Check 


When both transducer can be used as reversible transducers, 
the results of the reciprocity calibration can be checked promptly by 
exchanging the transducers' roles in the calibration measurements. 
Repeat measurement of the resonance parameters and transducer reposi- 
tioning is not required. 

Unless the two transducers have identical nonlinear responses, 
this second set of readings will confirm the accuracy of the assumption 
of reciprocity. It will not, however, measure the error in computing 


the reciprocity parameter J, which must be estimated independently. 


D. RESULTS OFTE ELİE ə 


1. Reciprocity Tests 


Tests were run on the three transducers in all pair combina- 
tions at the three accessible modes: n= 3, 4, and 5. Modes below n=3 
were inaccessible because the LC-10 transducers' transmitting response 
was too small for reliable measurements; modes above n-5 were attenu- 
ated into the ‘tails’ of the lower resonances, which obscured the ian 
power points needed to measure Q. Comparison ot a number of 
reciprocity calibration test results showed reproducible results to within 
0.8 decibels. However, the calibration results in Table IV disagree 
systematically with the values obtained for the same hydrophones by 
comparison with the USRD standard hydrophone (Table I!) by 1.5 deci- 
bels (Table V). 


2. Error Analysis 


The self-consistency error quoted in Table IV was determined 
by reversing the electrical roles of the two transducers without 
changing their positions; the difference in computed results for the two 
measurements is the degree to which the system is not reciprocal. The 
estimated error is the combination of the self-consistency error with the 


error in measuring the reciprocity parameter, listed in Table lll. 


36 


The primary source of self-consistency error is the difference 
in pressure amplitudes at the two maximüm locations used, which we 
believe is caused by attenuation. By careful selection of modes and 
locations, We were able to reduce this error (027 EGE ie ther lesser 
source of self-consistency error is the changing effect of the support 
tubing used for the transducers on the sound field as they are lowered 
into the water column. Insertion of the transducers and their supports 
into the water column displaces about 0.4 percent of the water column 
in volume, causing a similar change in the column depth and in the 
positions of maxima and minima. The change in the minima is the most 
significant error encountered in measuring the wavelength of the 
standing wave; because the cosine function is flat to first order, 
however, the effect on transducer positioning at maxima is negligibly 
small. 

Compliant materials associated with the signal cables and the 
transducers appear to degrade the reciprocal nature of the geometry. 
This error is reduced to very low values by selection of transducer 
location and by eliminating compliant materials (signal cables, for 
example) in contact with the water. 

Parameter measurement error is dominated by the apparent 
variation in the system quality factor over the length of the tube. An 
increase of 6-10 percent in the value of Q determined by the half-power 
frequency method was observed as the depth of the measurement was 
increased by one wavelength; our model for the compliant tube system 
does not account for such effects, and is therefore limited in the accu- 
racy with which it predicts the reciprocity parameter to about 0.8 dB. 
This is the largest recognized source of measurement error in the 
compliant tube reciprocity method. 

The systematically lower sensitivities measured by the reci- 
procity method relative to those measured by comparison with the H56 
standard are not yet understood. This discrepancy may rise from the 
inadeguacies of our model for the system energy balance, from the 
method used to measure Q, or even from technical difficulties with the 


comparison calibration. 
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3. Measurement of Wall Motion 


Since sound speed can be measured directly, it is not neces- 
sary to attempt the difficult measurement of the actual elastic strain in 
the tube wall as part of the calibration process. However, the degree 
of correspondence of the actual strain behavior to the model derived in 
the first chapter is an important indication of the validity of our 
method. 

The circumferential wall strain was measured by attaching 
constantan .250 inch strain gages [Ref. 16] to the exterior surface of 
the, tube at the depths of some of the pressure maxima located by the 
test transducer. A second lockin analyzer used as the sensing element 
of a resistive DC bridge measured the magnitude and phase of the 
strain oscillations sensed by the strain gage; by comparing the output 
of the strain gage with the magnitude and phase of the pressure meas- 
ured by the transducer at that depth, we computed the complex impe- 
dance (W,8) of the wall ring (equation(1.21)). The results in Table VI 
agree in general with the measured sound speed in the tube. Also, the 
indicated damping phase effects are reasonably small, justifying our 
approximation for the wall potential energy computation of eguation 
(57 
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IIl. UTILITY OF COMPLIANT TUBE CALIBRATION 


A. LIMITS OF COMPLIANT TUBE CALIBRATION 


The compliant tube reciprocity method can be extended to low 
frequencies by simply increasing the length of the water column 
enclosed in the PVC tube, or by decreasing the wall thickness to diam- 
eter ratio (equations (1.11) and (1.20)). The practical lower frequency 
limit is the length necessary to support the n=3 standing wave, which 


from equation (2.4) is 


N 
O 


1 mi CAD 


The upper frequency limit depends both on the onset of azimuthal 
modes, and on the increasing attenuation of the sound wave by the PYC 
wall. The former constraint breaks the plane wave symmetry when the 
frequency is greater than (0.586*L/R) times the fundamental longitu- 
dinal mode [Ref. 1/]. The attenuation problem causes excessive varia- 
tion of the measureable system quality factor with depth and reduces Q 
to the point where half power frequencies are obscured by adjacent 
resonances. 

Test transducer geometry is restricted to shapes that are small 
enough in cross section that reflections from their lower surface do not 
significantly affect the standing wave. The length of the reversible 
transducer, as discussed in the procedure description, should be as 
short as possible consistent with adequate transmitting response. (A 
flat disk transducer might be more effective than the cylindrical LC-10 
we used.) The test hydrophone can be much larger, if the transducer 
amplitude averaging correction is properly computed. 

Accuracy of the method is primarily limited by the effects of the 


wall attenuation. A derivation which accounted for these effects, and 
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for the lesser effects of disturbances to the standing wave caused by 
the test transducers and their support tubing, could conceivably 
achieve much lower calibration uncertainities. 

A significant advantage of the compliant tube method is the use of 
pressure antinodes as virtual rigid boundaries at which to locate the 
test transducers. Many of the difficulties associated with mounting a 
transducer in an real boundary which must have some compliance are 


thus avoided. 


B. POTENTIAL APPLICATIONS - 


Two applications for the compliant tube system on a practical scale 
have been suggested to the author: a field/small lab calibrator, and a 
large, very low frequency calibrator. 

The portability and flexibility of the test system used for 
compliant reciprocity recommends the method for general calibration use 
in small laboratories and in the field. The l-meter PVC tube provides 
calibration capability in the 1 kHz range; a 2-meter, lO-inch diameter 
tube would be more versatile and could reach frequencies of 500 Hz. 
The manipulator assembly need not be as complicated or massive as the 
one used in our experiment; any system which can position the trans- 
ducers within a few millimeters is adequate. Because the method 
involves resonance measurement methods as well as the techniques of 
reciprocity calibration, it may also be of interest for advanced instruc- 
tional laboratories. 

On a larger scale, a 10-meter, 30-inch diameter tube, which would 
fit upright in a two-story open bay, could reach frequencies below 100 
Hz and calibrate units similar in size to the Navy s TR-155 transducer. 
If current interest in very low frequency underwater acoustics indicates 
a need for an alternative to the traveling-wave tube system described 
in [Ref. 8], compliant tube reciprocity may provide the answer at a 
very low cost. 

There is also no reason we are aware of that the benefits of 


compliant waveguide usage should be limited to resonant calibration. 
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Construction of a compliant traveling wave tube calibrator with a non- 
reflecting termination or of a Helmholtz resonant calibrator using two 
compliant tubes should also be practical. The Helmholtz calibrator 
concept in particular has the potential of extending the range of accu- 
rate absolute calibration for the underwater acoustician to extremely low 


frequencies. 
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APPENDİX A 
GLOSSARY 


Subscripted variables in the text are represented using paren- 


theses, due to word processor limitations. 


A 
A(T) 
a 

EG) 

C 

D 
Ela 
Ed? 
EENG 


L(T) 
M 


D O О О 


Cross sectional area of the water column at rest. 

Active surface area of the transducer. 

Instantaneous cross sectional area of the water column. 
Parabolic curve fit coefficients. 

Sound speed in a compliant tube. 

Distensibility of the boundary wall. 

Transducer-driven, test hydrophone output voltage. 
Projector-driven, test hydrophone output voltage. 
Projector-driven, reversible transducer output voltage. 
Resonance frequency. 

The test hydrophone used only to receive during the reci- 
procity calibration. 

Thickness of the tube wall. 

Transmitter current in reversible transducer T. 
Reciprocity parameter, or acoustic transfer admittance, 
between two points. 

Adiabatic compressibility of water. 

Wave number at resonance = 21/A. 

Depth of the water column. 

Transducer s active length normal to the plane wave. 
Hydrophone voltage sensitivity. 

Number of quarter wavelengths below the water surface. 
Resonance mode index. 

Unit vector normal to surface element ds. 

Acoustic wave excess pressure. 

Pressure amplitude of the acoustic wave. 

Resonance quality factor. 


Inner radius of the tube. 
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R(m) Damping coefficient for tube wall vibrations. 


S Transmitting current response. 
S Surface area element. 
T The test transducer used both as a projector and as a 


hydrophone during the reciprocity calibration. 


U Volume velocity of the acoustic wave. 

u Particle velocity of the acoustic wave. 

V Volume element. 

W Complex mechanical impedance magnitude of the tube wall. 
x (0) Center coordinate of the data set used for curve fitting. 
x (m) Coordinate of extremum of the curve. 

Y Elastic coefficient for tube wall expansion. 

y (i) Sample set to be curve fit. 

ZTE) Transducer blocked electrical impedance. 

Z(m) Transducer open-circuit mechanical impedance. 

z Depth from surface of water column. 

z(0) Depth measurement offset. 

z(m) Depth of the mth pressure maximum/minimum. 

ፅ Sampling interval for curve fit sample set y(i). 

E Areal strain in the wall. See equation (1.19). 

£ Length averaging correction factor. 

8 Complex mechanical impedance phase angle of the tube wall. 
À Wavelength of the standing wave. 

p Instantaneous local water density. 

p Ambient water density. 

$ Transformation factor. 

x Sum of squared errors in a curve fit to data. 

ሠ Angular frequency at resonance. 

č Areal expansion coordinate. 
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APPENDIX B 
PARABOLIC SMOOTHING ALGORITHM 


The algebra of deriving formulas for the best choice of coefficients of a 


nonlinear equation, such that the sum of squared errors described by 


12 
=D = E İytə - il (B.B 


is minimized , is described by Bevington in [Ref. 13]. The minimum 
for the function defined by equation (B.1) is defined by the selection 


of coefficients which solve 


(5:28 


where C is any of the coefficients of the curve to be fit to the data. 
The simplest function that will predict maxima and minima is the 


quadratic equation 


EE E 


(B° 
The single extremum of this curve occurs at 
X m — Xo = -C,/2C, (B.4) 
with respect to the coordinate origin located at x. Let the sample set 


be seven measurements of the function y sampled at regular intervals 6 
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Y1= y(Xo +10) ; i=-3,-2,-1,0, +1, +2, +3. (B.5) 


Evaluating the best fit function's partial derivatives with respect to the 
three parameters for such a sampled data set gives equations (B.6) 
thru (8.11); summing and cancelling yield equations (8.12), (8.13), 
and (B.14). These analytic results can then be used with equations 
(B.3) and (B.4) to obtain the location of the maximum and its ampli- 


tude, or the smoothed value for the amplitude at any location. 


2 (B.6) 
X-3 = 96 C2— 36C,+ Co- Ve 
X -2 = 46°C2- 26C,+ Cor Y-2 KE. 7) 
X-1 — 6°C2- 0121 Camy- (B.8) 
Xo= Co- yo 
(B.9) 
X +1 = ó^ C ÓC;+ Coan Y +1 
2 (8.10) 
X+2 = 46°C,+ 20C,+ Co- 75: 
y+s = 96°C2+ 30C ¡+ Co— Y+ (Bol) 
2 (512) 
840 C2 = 5(y.;+ Y -3) - 3(y. y-1— Ay 9 ) 
286, = 3(У.з— у-з) + 2(ys2— Y-2) + (yar yi) p 
2100 = 7Yo ES e6(y. + y T SY +2 +y-2) = 2(y+z+ Y -3) (B.14) 
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APPENDIX G 
ACTRE: HP-85 BASIC PROGRAM FOR RECIPROCE EE 0 


A. WHAT THE PROGRAM DOES 


The Automated, Compliant Tube Reciprocity Experiment (ACTRE) 
program provides the operator with a menu of soft key controlled 
"tools" for resonance parameter and calibration measUremeniuse [ts 
design was based on the need for maximum flexibility, rather than 
maximum unattended operation. 

All of the measurement "tools" are subroutines called from the 
HP-85 function keys by the operator. The major procedures are: the 
peak frequency location algorithm (KEY 1); the spatial optimum location 
algorithm, which also computes wavelength and sound speed (KEY 2); 
the Q measurement algorithm (KEY 5); and the reciprocity calibration 
measurement sequence (KEY 7). Supporting procedures provide the 
actual control statements for remote switching, function generator 
setup, manipulator positioning of the hydrophones, receiver output 
measurement, and graphics output. Any procedure may be terminated 
by calling the process clear sequence (KEY 8). Comments on the 
detailed structure of the program are includéd with the source code. 

The program was developed to fit the demands of the experimental 
apparatus,the experimenter, and the research goals of this work as well 
as the limits of the host computer permitted. Time and core memory 
constraints severely limited efforts to make the program friendly or 


accessible to the uninformed. 


B. TES? EQUIPMENT BUS CONNECTIONS 


The code is designed to use the IEEE-488 bus-controlled test 
equipment described below. All of the primary measurements are 
made: with bus-controlled equipment, except for the hydrophone output 
voltage measurement by the EG&G 09204 lockin analyzer, whose output 


signal must be converted using the HP 3421 D/A unit. 
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Bus assignment numbers used in the program: 


E 


705 
706 


707 
709 


714 
120 
[ru 


HP7470 hardcopy plotter. 
HP59307A remote controlled switch: 
A switch: reversible transducer; 
B switch: receiver input. 
HP 3314 function generator. 
HP 3421 Data Acquisition/Control Unit, with Opt 020 
board in slot O.. Channel assignments: 
O: Manipulator motor drive up actuator 
1: Manipulator motor drive down actuator 
9: 9204 lockin analyzer magnitude signal. 
HP 59307A remote controlled switch (transmit side). 
HP 5316 Frequency counter, in count totalize mode. 
HP 3478 digital multimeter, configured as an AC volt- 


meter. 


VARIABLE ASSIGNMENTS 


| 


FNH(X1) 


FNK(KO) 


FNL(X0) 
FNM(XO) 


FNR(X0) 


Defined Functions 


Length averaging correction function (equation 
(ESSE 

Cycles KO through the indices 1,2,3. Used in plot 
scaling algorithm. 

Returns decibel value of X0 = 20 log(X0). 

Returns voltage sensitivity in dB re 1V/uPa (equa- 
EE oy). 

Returns value of XÜ rounded to three significant 
places; reguired for compatibility with HP3314 function 


generator amplitude command. 
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2 


Variables 

A (7) Amplitude sampling array, used for curve fitting algo- 
rithms. 

A2 Amplitude ratio of hydrophone output voltage to 
transmitter input voltage. Measured by the Amplitude 
Read procedure. 

A4 Half power amplitude in O measurement procedure. 

AS Lockin analyzer output signal to D/A unit (709) 
channel 5; ranges from 0-1 V to indicate O-full scale 
magnitude reading. 

A/ Resonance peak amplitude, computed from curve fit. 

A8 Preceding measurement  storage--in О measurement 
procedure. 

BIG) Plot scaling parameter array 

B1-B9 Plot parameters 

CES Sums used in wavelength measurement least-squares 
procedure [Ref. 13]. 

CI ur ona 
C2. Sumo NOR 
Co SEmMEor ZE 
C+ Sumo 65 
C3: Sum OT sm 

C6 The quantity A in the least-squares fit. 

GT Resonance sound speed (m/s). 

CS Error estimate for C7 (m/s). 

D1 Water column depth(mm) 

biz Water column cross section (mm?) 

D3 Water column volume(m?) 

DS Length of T (mm. 

D6 Length of H (mm). 
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E 
E2 
E3 


E(PT): Output voltage of T , P driven (V). 
E(PH): Output voltage of H , P driven (V). 
E b 000 olaca 66 FI. PF driven (V). 


All frequencies are in hertz(Hz). 


FO 
Fl 
F2 
F3 
F4 
F5 
F6 


F7 
F8 
F9 


G2 
GS 
G/ 
G8 


Center frequency of sample set. 

First peak detection / Lower 1/2 power point. 

Second peak detection / Upper 1/2 power point. 

Peak frequency estimate. 

Lower search boundary. 

Upper search boundary. 

Drive frequency; argument passed to Amplitude Read 
Setup procedure. 

Resonance peak frequency. 

Error estimate for F/. 


Stepsize. 


Receiver preamplifier gain. 

Lockin analyzer output expand setting. 

Calibration attenuator value.” 

Drive boost factor for preamp calibration test in 


Calibration routine. 


Length averaging correction factor for T. 


Length averaging correction factor for H. 


əl 


18 
19 


JS 
ni 


KO 
K3 


ES 


M1 
M2 


NO-N2 
N9 


Q7 


RO 


General use 

Size of plot files XQ), YQ. 

Selection index for Operator Input procedure. 

Number of time constants TS to delay before measure- 
ments. 

Select index for Switching procedure. 


Wavelength measurements count. 


Reversible transducer(T) transmit current (A). 
Reciprocity factor (V*A/Pa?). 


Water density  (kg/m?). 


Plot scaling array index. 
Pressure maximum/minimum position (m). 


Hydrophone voltage sensitivity of T (V/Pa). 
Hydrophone voltage sensitivity of H (V/Pa). 


Position counts; used to compute manipulator motion. 
Number of 1/4 wavelengths from surface corresponding 
to 


Measured quality factor of resonance. 


Position encoder count ratio for down "motis 
(counts/mm). 

Manipulator drive rate (counts/sec), down motion. 
Manipulator drive rate (mm/sec), down motion. 


Allowed projector drive adjust error. 


R4 Allowed peak frequency detection range. 


R5 Allowed min/max detection range. 
R7 Q measurement coarse stepsize factor. 
R8 Manipulator up/down motion rate ratio. 


The S variables are all control flags/indices. 


SO General use. 

S1 Read enable flag. 

52 Manipulator motion stop flag. 
S3 Plot option enable flag. 

S4 Fast display refresh flaş. 
5. Drive voltage select index: 


0: use J11 projector drive voltage V4; 


1: use piezoceramic drive voltage V3. 


50 Automatic drive voltage control enable flag. 
— Resonance peak detection index. 
5 Plot execute flag and abscissa select index. 


All times are in milliseconds. 


T? Manipulator motion time adjust. 

TS Lockin analyzer time constant. 

HT Additional delay time before measurement. 

[6 Size of additional delay T/ to be used when required. 
UO Curve fit amplitude at center of sample set. 

UT Curve fit slope*stepsize. 

[12 Curve ft curvature *(stepsize)*. 

U9 Stepsize used in last curve fit ( for plot procedure). 
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All voltages are in volts(V)rms, except where indicated. 


V1 
V2 
V3 
V4 
V5 


V6 
77 
Va 


WO-W3 
X (50) 
Y (50) 
XO 


Transmitter input voltage. 

Hydrophone output voltage. 

Piezoceramic drive voltage(Vp-p). 

J11 projector drive voltage(Vp-p). 

Lockin analyzer fullscale setting; input as milli- 
volts,and converted to volts. 

Function generator drive voltage (Vp-p). 

Transmitter input voltage at resonance. 

Curve fit amplitude at F7, returned by 5-point 


Voltage Measurement procedure. 


General use working variables. 
Abscissa plot storage array (Hz or mm). 
Ordinate plot storage array (mV rms). 


Reference x-coordinate for curve fit plot (Hz or mm). 


All variables beginning with "Z" are in millimeters(mm). 


210 
ZO 
22 
Z4 


Z8 
29 


Sampling positions for A Measurement procedure. 


Manipulator table position. 


"Estimated position of pressure minimum/maximum. 


Position of center of sample data set/ operator-ordered 
position to move to. 
Measured stepsize in À Measurement procedure. 


Ordered stepsize. 


3. Character String Variables. 


AS(128) 
MS(20) 
ssis) 

WS (32) 


Title; used to identify printout and plots. 
Code identifier. 
Remote switch position identifier. 


General character input. 
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D. SOURCE CODE 


1000 
1010 
1020 
1050 
1040 
1050 
1060 
1070 
1080 
1070 
1100 
LD 
1120 
LESO 
1140 
18-29 
1160 
1170 
1120 
77: 
500 
110 
m 
1250 
1240 
PESO 
10 
eg 
5. 
od 
1 300 
5:7 
70 
1220 
5. 
20 
75: 
55: 
1: 
400 
1410 
420 
1430 
1440 


1 


EEN aetres/ 

) 1070136860) Lonpiyint Tübe Reciprocity Experiməsnt 
! M.B. Johnson, Dec 84 

| 

፡ፆ፡ ከ ከ 01: 0 0111107 000 18 DES 

A ARCSIN Sf05] 
A$="ACTRE-7 mod 8 UE CERA DISFRE 
ringe ons 

ERT DE 

Eer = Ki MOD 34! 

ксри ДИ ју 1:041/111) 


Deb GE Meer dn 

EE EE EE +igures for HP.:14 
ири 

ከመ 1 ከ ማንን 

ከሴጨ፡ን65ውሱኙፕ 7 TREN ስ ለ1ን))ህሀ E Wi=W1 1000 6 5010 1140 

IRA ህያፐ ኛውን ርክብ (1710 € Viİ-V1:10 € GOTO 1140 

! end loop 

PNR Hu 19 1007.) 

FN END 

"phot “ünetlons 

wevear Pane initvalrzation 

2517177075 EE, E (5. C/7,La,Li,. 
REAL D3,E1, 57576 et ER. 
REAL NO, is ,N=, UC, Ut,U2, V1, 5. 25077 uM E „205 
SINT (3 O pq ህብ 0 ትጋ FB,F7.65,37,39 


F5, ’ ’ ነ | 
“Eu KE дә 577: 12, T5, EINE VAM VO. n 
MOSM AD 1030) 
| NRPOERSKE 27, ከ. E. 
5175::50721 (82.92.34, 2 
E EE 2 
205 obs! 
Sup иҹ̧) 
7:4) 9271021700 በበ... 
7177: 7:5 5:67 8 5600 
:77FP: 5775 ):0077 52.0 RO G aS=] © 121-198 © move seals 
ти есас 1 узЕ error 
ከን L OL በ0238 ә arror 
& GCSUB 20850 
EE 060 * posicion ን 
Sst Geese eden th. tenperatu 
ከ ES at) 
(2122 Aa poe scalrng parameters 
Т8=0 ! read delay factors 
ንዐ ወ በ ህከ water density 


кл 17 
e 


7 


re 40 –4Ф es 
tt i " 
t9 O) — 


газ 


ii 
AJ CH -— +— 
13 
O 


p-> 
tu 
(a 


C m ~nd 2 = 2 


–-- 11 23 


ә 
i 
~O 
(o газ 
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115157 begin: equipment 581 ህወ====2252222222፳2222222መመ2=ሙ2222222=== 
1460 OUTPUT 709 4"RS' d clear bK E 

1470 QUTPUT 707 ¡"PR" 1? clear 3211 2 1) ии 

1430 OUTPUT 725 "1277701 ህር. ወ 1, 

1490 QUTPUT 7107 s E DOS 

1300 QUTPUT 709 s"FIRORAOZINSLSS; Ti” i| —a7d recs ESQ 

1510 OUTPUT 720  "ЕМ12" TLA MEIS 

1520 QUT? ON 

1350 1391 0 GOSUB EE 

1540 1 end: equipment Sebtup=suss==s=zmananaszassaszszmnssnannanasanasunmiuas 
1530 ! begin: fresh slate for Processing ===========а==лиҹ=тата--- 
1550 13,19,N9,N420 ! plot file &w1/4 file sizes ,w1/4 count 

12970 C1,02,635,09, 05, 0070 0 WAG) ENENNNESUME 

1580 =7,F8,C07,C8,07,V4=0 

1590 Hİ 000000 0 5631 

1500 V4=,0001 6 V3=10 8G T7259 

1610 Söe1 ! ADC on 

1520 IES TOR 50300 0100 5: 

1630 GOSUB 1800 ! set main menu 

1540 ! end: fresh slate for proceszingessususmuuxzuzuzsummuzumusn 

1550m beginipreczssinę ann x. .....ı...-..-.-.............5- 

1660 ! hold here for operator menu selection 

1570 CEEFIR 

1439 DISP USING 8590 ; AŻ,F7,F8,C7,C8,07,FNM(MI) „FNMCM2) 


1690 IF S3 THEN DISP "plot on  ";ELSE DISP "plot cfi —"; 
1700 IF 85 THEN DISP "ADC on" ELSE DISP "ADC off" 
1710 DISP S£;" 10=":20 


| ህር... 
7300508 101107” 


1750 IF S^ THEN 1740 ! fast refresh on return trom ev 1207700 

tee 

1760 543120 

L770 G07 ole 50 

1700171 Drocessgina 251395።፡፡፡፡ 77፡2 25:22 522 2:5:፡ 2 5:2:512%.. as 
` 
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1800 ! begin: MAS === === =2532= 230 Oo a SIS 
7170116757 0, En” GOSUB 2650 
1820 ON KEY# 2, CHRS ( 11) GOSUB 4150 
1830 ON KEY# 3,"move" GOSUB 7990 
ee ON REM, oor’ GOSUB 2040 
1550 ON KEY# 3,°Q" GGEUB 3390 

ር ከ ር በ isa pe: bDSUDp 6120 
KERALA GOSUB 5170 
5757573117 clear 5010 1950 
1890 S4=1 ! fast refresh 

1900 RETURN 


፡:.፻0 ፲ 20:ሀሠ::ሀ291[ሀ. eN EE Soo anno 
po begini e ear moate menú 
1730 ! provides exit/retresh ки 


Ee erte reiege DEE TIMERE Z 

2777711 0117707 7 DP NO? E OUTPUT 709 ("ОРМЫ1" 
5255::1 | Flag reset 

TO 255! plot dumo 

AE main" 6DSUB 1800 

ከ ንውኤሌሒ) ሲከ Lu “109” 5070 1550 

21 7 ON 5671770 quit" GOTO 8640 

PRINT ED CLEAR/RESET” 


7C:(VU7) )) si re”resh 
7 7)::5211 1000 )İ entar processing stage 
eno eC Ee MENU m O m a a o”: —- 
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2050 ! ======PROCEDURES======== 232 SR SEE EEE PV— EEE ии ии 
2060 ! begin: operator input procedure (HP-HUMAN bus)]=========== 
2070 Ir 1450 TEEN 2059 ELSE ey 

2080 ! menu display 

2090 CLEAR 8 DISP USING 8500 € 1155 11000 ии 

2100 INPUT WO 

2110 IF WO?9 THEN ZO=WO 6 WO=0 ELSE 14:40 

2120 IF 14=0 THEN 2630 | return te naim mena 

2130 ON I4 GOTO 2140,2200,22850,2310, 2580 2-200 SIL LLL 
2140 ! 3204 LIA settings(i4=1) 

2150 DISP "LIA fullscalelmV),time constant (ms ero ду 

2160 IİNPUT Vo,T3 52:02“ ከ ከ ከሮ 

2170 DISP "delay: 8 time consts,extra timge(ms)' cd 1872011 ии 
2180 PRİNT USING 8510 7 1000270000 ME 

2190 601175 

2200 | water depth,templi4=2) 

2219 DISP "column depthimm), temperatura L7";Y NEED WARAN 

Seo ny. 

«2900 “mel As DNS GD ANNĘ 

240 PARIS O Oo DR) 

E EE 

240 ! drive voltages(i4=5) 

2-70 DIS “drive for ኮሦ/1| УРЛ ви ии. 

2280 IF LEN(WE)>0 THEN V4=VAL (WH) ü VOSVAL  Tİ”UBUU ди 
2290 PRINT “P ürive:":V3:"VU" € PRIM o o EN O O D 
2200350006 

2310 ' comments (i424) 

2.0 0165 "comments" 

2330 INFUT Ni 

2240 IF LEN(43)20 THEN PRINT "BE". © 000 2060777” 
2070705) 

-s50 l API GUS Constanze si) tao) 

2.70 CLEAR € CISF "enter new values,then CONTINUE" € FAUSE 
ес 


ርያ 50761 2500) 


2400 
2410 
2420 
2430 
2440 
2450 
2460 
2470 
2480 
2490 
2500 
2210 
28.20 
29350 
2540 
Dos 
2360 
2340) 
EA 
msc) 
2600 
2520 
250 
2550 
2540 


120507 426) 

AF=AFT1,18] 

DISP "date"; INPUT W$G A$-A$&$4$ @ GOSUB 5070 

DISP "run #" @ INPUT W$e AF=A$LW$S 6 GOSUB 8070 

! transducer chars only(i4=9) 

A+=4$[1,64] 

H$z"T is:" @ DISP H$;@ INPUT W$8 A$S=A$8GHSŁNS A GOSUB 8070 
2515570006 6 0120 mm) Capa cancee ab "6 INPUT 0S,Y3 
H$z"H is:" 8G DISP H$;8 INPUT W$8 AFSASAHSANS 

DISP. "H active lengthimm) ';G8 INPUT 05 

PRINISUSDNG. 9520.5: AS,05,Y3,D6 

ie) 0 10007 00 € 0206 1000 ! SI conversions 
SOTO SD 

' plot enable/disable(i4=7) 

232eNÜT 52 8 [3:0 

DEC ሎክ От ЕЕБЕ 2650 

2726 15 705 

7757:5:577070 007 E 22,122,27,7/ Y PEN A 
5070915) 

! position (i4=8) 

Mi pesicion= 70:8 INFUT Wf 

İF LEN(WS)>0 THEN ZO=VAL (UF) 

20100 2350 

14:0 € S4=1 @ RETURN 

' end: operator input proczduresssntumacznuzmuzzmacz zu 


E 


2650 
2660 
2570 
2680 
2690 
2700 
2710 
2229 
275: 
2770 
272) 
2750 
2770 
2780 
27: 
2300 
2810 
2820 
2830 
2840 
2650 
2850 
2070 
2880 
2890 
2900 
2510 
2000 
27320 
2940 
= 1 
2950 
90.0 
2983 
2790 
3020 


! Begin: peak frequency məasurement procedure=s=======z=-su=s= 
PRINT USING B4 509 ии из 

CLEAR € DISP “peak frequency measurement" 

DISP "estimate,search range,stepsize";8 INPUT F7,NO,F9 
F4=IP(F7-W0) 

ት ር ርህ ከ111 

! documentation 

PRINT USTNE 847073 ጊጊ... ክን 

! warmup reading 

IF SS THEN VES VE 6 66” .) 

FA=F7 ( T7=T8 € GOSUB 6990 JEE HESE 

DÜSUB Ee 

IF 53 THEN 59=1 € 60SUB 3770 | Gegln 17:10 TA E 

| begin tinitialize Samplinopas 

FEOS PESE 0600077) 

A d sa 

GOSUB 7210 ! read amplitude A2 

IE 1207 THEN FózFó+F9 © GOSUB 6990 ! xmtr on next 0440777 
Airey ene 

NEST AT 

GE a ER ER e E 

' begin:resonance search 

SI, F1,F2,F7=0 | peak detect tlag and de:5660000 7 

FOR poses 1077000 Oo 

| +0 is the freg at the center ofthc samp vec. 
Fa£=sF0+4+F9 € 52505 6970 ! «mtr en nest Falacia 
ርህ ከ ከ ር at. 7c 

541141. 41.72.4109. 51211621. 141.70 

DISP USING B380 „ ONE 

KEE ሄክ: 

EE 
Ui=(3é(At7)-AL1)1)+2 lA OA JETI 
| peax frequency 

505057 

Ri E A C EES 

IF АВ50Р9-20)6294*2Ф THEN 000 07 


i2 5 2 
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3010 ! peak detected at center of sample array 

3020 ! peak amplitude 

2177 Ao 0 060 (A (21 +A(6hh-2x(A(7+A(1)))/21 
3040 87zU0-U1*U1/(4*U2) 

2050 S/=S7 tle V/=V1 

S060 iF 55 HEN VSsVS ELSE V4=V6 

GO EN RS ELSE F2=FS 

7751110660 0007 o 00 Vİ, İ0Ü00sU2, 1000:V1 
3090 U9=F9 € XO=FÓ0 

£ 00 5101702177 | 

22110 eqndsaesk detected Block 

3120 ! begin:peak not detected block 

3130 S720 ! unsst peak detect flag 

2140 ' end: peak not detected 

2100 1 end:tests at FU 

310 DTS. 

27?) MEN Sta0 ELSE S240 | end search on two detects 
¿3186 | gush down sampling array 
571577) TO ë 

“200 ACIO)=A(IO+1) 

510 NEXT ID 

50 50202 Tele 6 A(7)=AŻ 

rule (1 FO 

57200 ' end search 

5751175. 768 GQZUB 3990 ' show plo! 
1257117557.) LU?” 6070. 32/70 ELSE 2310 
2270. IF F1*0 THEN PRINT "search failed" @ GOTO 3350 
7289 ! one gcod reading 

፡ ን ተ! FSF3 

77 001020 

57777 ewe good reudings 


| ለ” መጫ*ፆ .,. 
» P > "ei A? 
= = ሎፍ= '— am 'ፃ 


DI 


12757) e average ol mO good ra., dings 
2430 | final results 

፣ ከ UslNs 9310 ; F7,F3 

5505262“) * fast refresh 

E530 RETURN 

2370 ! end:center freg ant 07ያ0899የ2=፡5፡2፡=:፡:ዊ፳ዊ2፡22255ች:52=2፡2:፦::2 
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3380 ! begin: Q measurement procedurasz====c=====s====sx=zzzuuusr 
3390 CLEAR @ DISP "@ measurement" 
5400 DISP "center peg ee 
3410 0720 
3420 PRINT USING 8570 3 TIME 20 
3430 GOSUB 4950 ! measure peak voltage v8 at f7 
3440 A7=V9/V1 
3450 ! set search limits 
3420 ትተ. ህህ AA 
$470 F4=MAX(0,F77204F7! 
2480 F5=F7+20+F9 | 
3470 DISP "sneak (0.10005 25 
3900 A4z87*.7071057811825 1 1/2 2012770107 
2510 DISP “hal? poweri 20 DEM 
520 ! Baginisecarei for lower I TEn 
(000455 MR E 
2940 58247 
3330 Bola aloeg to tinaoa 7 
60 Et:kbo-ev 
2573 IF F69F4 THEN GOSUB 5990 ELSE 4090 1 /m0 7 eee 
+ 0 20505 72:16 |! measure AŻ 
390 E ebe RE 
2400 ' ChEGE Wer ‘coun ከከ ከከ... 
ool) IF 42209 MEN (62) 6. ን... 
2620 ' Uphill, lew 1/2 5600560027 
3530 PRINT USING BE, EEN: 
7640 GOTO 4110 
2600 İ downhrllersset checks 
2660 AE=AZ 
SETO LIF Eege THEN’ S500 ЈИ 
2250071 end een 
2370 ! begin losn: pracko NE n 
2700 1Е РЫҜР7 TREN FósFó+1 4 G0305 36077. 
3710 AB=AZ ! save last sample 
3720 GOSUB 7310 ! new sample 32 
2730 DISP "aq'":152-A4)*1000*U1 
“26078 A2>=44 THEN 3/:075L55 ብ... 
21309 | end loop. bracheted “1112727797 
A A A A ey ie ae. 
270) RINI 1016016307 s 


N | ` ኣመ ” = 
eee an “6072707626 || 
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770 
3800 
3810 
370 
3850 
3840 
3830 
3940 
2870 
3980 
98 71 
3900 
52:10 
75 
27 9 
3940 
70) 
577 
2770 
2980 
2990 
4000 
4010 
4020 
4020 
4040 
4050 
4041) 
4070 
4030 
4090 
4100 
21:10 
1120 
4120 
4140 


Eno le power pt F2 


FĘ=IP(F7) 
48:47 | check value 


een la fóta 


Fó=F6+F9 

IF Fa€F5 THEN GOSUB 6990 ELSE 4090 
GOSUB 7310 !| measure AZ 

MS 42700) 000071 


check downhill motion 


15:52 RE 0) 3720 


! 


conil 1/2 coscured 


PRINT USING $520 ; F4,82,A8 @ COPY 


GOTO 4110 

! downhill-rëeset check val 
АЗ=А2 

IF A2<=44 THEN 3950 ELSE 5820 


i 
* 


end loop 
"Eon oop. bracket 1/2 point 


12755 AEN Fó=Fó-] © GOSUB 6990 ELSE 4090 
AB=A2 ! save last sample 

БИР 7510 1 new sample a2 
HEEN ` 

IF A2>=A4 THEN 4020 ELSE 3960 
ብር I085.bricksted-interpolate 
F2=F6+(A2-A4)/ (42-43) — 

ር ከከከ ወከ 110430 010 e 

r= 51)” quğliity Factor 
S4=1 ! fast refresh 

"..ከ ከብ. ፣ IND 8707 

34=1 1 fast refresh 

RETURN 


begin:error processing 


ANS NG BDO + Fo 

Eu 

GOTO 4080 

' end:@ measurement ри осефуласиз=алтлилла: 
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4150 
4160 
4170 
4180 
4190 
4200 
4210 
4220 
4230 
4240 
4220 
4250 
4270 
4290 
4790 
4500 
4510 
4320 
455 
240 
4250 
4350 
7: 
4580 
4190 
4402 
4410 
4420 
4522 
4440 
4450 
4450 
4470 
449) 
4499 
4300 
4510 
4520 


4550 


50 


1340 


4350 


! begin: wavelength measurement proceduresszszzzzzzrz-zzzzz- 
CLEAR € DIS? "minimas posu tc 

IF SS THEN V&sV3 ELSE VA=V4 ! select T or P drive 

ЕбзР7 8 Т7=Т8 € COSUB 5990 | xmtr eno 
GOSUB 7310 ' warmup reas 

IF 55 THEN 5752 @ GOSUB 3770 | phot sto NEM 
DISP. "positions 20 

DIS? "stepsize በ ከ ከ ከ16. 

PRINT USING 3727771157 

| Bagin:ial 218112 205 4A6 a 

ЕОК 12 907 

GOSUB 7310 Һәзи әтри иди. 

24-10 6 GOSUB 7510 ! move down 29 mm 

2(121*224 8 A(I2) 252 


(4 


NEUES 

ип initia li Erana 
S7=0 

' begin: search loop 
14=1(4) ! array center 


1 uma 0017507710 


U= (ILALA -IAA Ia G 1) 
Ule(:(RA(/71-5(1116 004 eam Ge 
UO=(7FARIL LES CAC) TI EE (7) AiL: yopa 


19=18 8 *0524 4 A7=UN-U1+U1/(4+U2) 
12=14-18411/(222) ! predicted optimum 

DISP USINE (85300 24 E 

1567 OR aBS(22-243)4RZs:9 THEN 4470 ELSE 4460 
| hit- document 

1127 THEN 4280 l የ የሚ ው ሜስ 

S7=1 ! set search end 

Zi=Z2 ! save peak estimate 

' continue search 

' push down sample array 

FOR I0=1 TO £ 

116) =2 (10-11 € “1101212017 

NEXT TO 

BOSUB 7310 Pro E) 

2170520 8 A(7)=AC 

IF 87 THEN PRINT USING 8720 y 22,12900x9/*V1 UI GS U > 


^nT KGR 
8010001 D 


' end: search l66Q 
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4560 
4570 
4530 
4590 
4500 
4510 
4520 
4630 
4440 
4630 
4560 
5670 
1680 
1690 
4700 
4710 
17: 
1720 
3740 
ጳ7፳0 
4760 
4770 
1722 
6790 
4220 
4810 
422; 
4520 


AZ 
“77 


<Q O14 CO 
~E F} O LA 
F ES GO A 


. 4 


¿890 
4950 
4910 
4520 
17.30 
4740 


да -— > -> Za 
oo 03 


LO 


| print final estimate 


PRINT USING 8730 ; 22,1000*A7*V1,U1,U2 


12=(22+13)/2 ! average 
PRINT "average=";22;"mm" 


EE EE EE show plot 


! begin data reviaw 
IF S& THEN 4630 ELSE 4640 
1557 JREN VSSYS ELSE Y4=VG 


store 


ADC results 


DISP CHR$(11):"/4 count=";N9+2:8 INPUT W$ 


KE AO OENGO LE 7-7 6 


HÜZVAL (45) 


TF WO>O THEN N9-W0 ELSE 4910 


' beginiaccept measurement 
55: 

ዘ 5271000 | meters 

' begin: sums 


55507 “ር 

5565001217 7” 

2::77 57: 7 

65001 | «xy 

ee 07 уб 
5I9*C2-DL DI ! del factor 


-> የ3 


ond: sums 

' begin:ccapute sound spd 
ZRYW JIEN “610 EL9t 4879 
Meie cr more 
ızı 1:00 Lo /660 ! 5 
"5545 00/15 1 
77 P5:757:7 

1541772 THEM “060 ELEEMƏYO 


Ye 
* 20 
mmo 


С3-=44*5028112521244)04443Ҹ1441222-2 
1$:9/C65Y *F7 


| document 


2 
rcegt 


ግ 


NOR ge 


PRINT USING 369% ; 19,07,09,1600€1 


COTO 4930 

! end : accept measuremt 

| Begin: discard measuremt 
PRINT "mmt discarded" 
27570 RETUBM 


' @nd: wav2length measurement 
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6670 4680 


መጻ ~ 


^ ' . E Oe : 
ASMHİ 222 un. š 


ат nem mg x Er wh ah na w ma pz w ua <. መ ጭጩፎ ‹፡“ 


4930 
4950 
4970 
4980 
4990 
S000 
5010 
5020 
3050 
3040 
“də 
3060 
3070 
5080 
2077 
3.100 
2110 
SREB 
2: 
5140 
130 
3160 
3170 
5180 


! begint -pt peak voltage —X measurement procedures ን |. 


IF $5 THEN V6=V3 ELSE V5=V4 ! select T or P drive 
F=F7 6 Т7=Т8 € GOSUB 45990 ! xmtr on ,long delay 
GOSUB 7310 ! warmup read 

| begin:five-point mmt of peak amplitude v2 
50=56 8 55=0 ! ADC of? 

F5=F7-2 8 GOSUB 4990 ' set xmtr Fn-2 

GOSUB 7310 "pead v2(Fn-2) j R Dp 

FósF4+1 @ GOSUB 5990 ! xmtr on Fn-1 
Vaz-(3:02700) Flies cern 

SOSUE ብከ. ከ ከ ን 

FósF4+1 ፳ 50508 6990 1 xmtr on Fn 
VB=VB+12*V2/35 ! 2d term 

GOSUB 7310 ! read at Fn 

FózsFó+l G GÜSUB 5990 ^ xmtr on Fmi 
V8:V8*17*V2/53S ! center tera 

G0SUB 77210 !' read at Fn+l 

Fóázró*1 à GOSUB 6990 ! xntr on Fnt2 
V8=V8+12+V2/55 ! 4th term 

505409 7909 Pİ peda st Pn 

V8EZVB-2*V2/25 ! result at f7 

96:90 ' restore ADU 

RETURN 


! end: peak voltage ӧфӱӧсла=лазыневнзвтаниад олн ии) C 
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! begin: calibrate hydrophone proceduressszzzszzzzzz2225£50727 
412C7/F7 ! wavelength 

Hİ-FNH(Pİ:03/V1) ! length average corrn for T 
Н2=РМИ(Р1206/41) ! length average corrn tor ዘ 
ANAIS IDAS ከህ ከ ን NL (E), Fr NL (Hz) 

| beginspreamp calibration 

le EE pe 

18-4 € GOSUB 8120 ! switch 

DISP “preamp calibrate" © WAIT T3*J 

2501001 o 6000702 76 1 “13 8 60565 5990 : xmtr on 
GOSUB 7310 ! read LIA out 

G2=G7*A2 

PRINT USING $740 7 FNL (G2)  1000¢4V1 

' end: preamp calibration 

' Tid measurements 

155.78 50508 8120 ' switch 

Bakes: 6 INPUT W3 

DOL UET gpt mat ef ve 

ከ ከ G Tl Eth 

JZEVIEŻRFIXFZEŁY, | amperes 

ት ከሁዜፕ ፕ የለይ፡ቻልያታተታ ጅእከኪ ስለ U1,1000475,1000000ZE5 
A T; 

ከ ከ ከብ A 

Elek 

EK NA 

5:7: (750705: 

3450 Fó6=F7 6 VasV4 8 T7=T8 8 60508 6990 i «mtr on 
2:500 5USUB 7210 | warmug reading 

2470 ! begin: adjust loop 

EE 1 xmtr en 


CO OOO O Ch ጋ © © 


(ባ £n ርባ (ባዳ cn £n on CI CA ON (A 


“3 O 4.> CI фә LA r) ** O -0 DI ሩ-4 ሁ (ባዳ ፦ ፈ-4 ዮጋ = CH –-Ә 
© ҸӘ © © OO © C» €» © 


o 


(n cn ca cn can CA CH LI LN CR cn 
ዝብ» M od b Ll i Ci £4 GA CA LĄ 24 CA CA rJ 723 r5) £2 PF) FI PY FO ያ.3 ፻›3 ፦፦ 


1 
) 


e 
m 


1 


са сл са En 
KC 


ES 
-P 
C» 


Ü UZ Z101:: ርር ን 00 0 ' read v2, correct ror nream 
ማንስ 7101006 2710? 1000%V6, LODO V2/62, 10000000 7 
5:77 16-00 060 0 HE Vo: NN) (6 £3/£2) @ 5219 5470 
2/0 7 eng adjust lech 

7.0 /4=/6 

7”O 1.11” (rTives"?/U4:"V" 

Sza". end FP drive adjust 

2777070210 GÜ | ህሪ mat 

s ክህ የ 2/5032 1) Eph 

2777710107 cu”) (000000x22:"uV” 

3077) ena PON 
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5600 ! begin P3T 

3010 DISP "Р)1" 6 185-1 š 6051000000” 
3620 GOSUB 4930 ! vB mmt 

5650 ከ1መንህ 1. ከከ 

5640 PRINT " Ept s";1000000*E1; "uV" 

5550 ! end P)T ants 

3660 ! begin sensitivity calc 

96/0 J7=PTRDSSE7/ ROSE 77800 

5480 PRINT USING 8750 ; F7,07,07 

39690 FPRINTI2USI O 77 UE 

3700 M1=SOR(J7AE1+E3/(E2*J3)) 

2710 ሸሬ=መመ ARN M 

5720 PRINT USING 8780m EENM 11 E 

3/30 ! end sensitivity comp 

3740 $421 8 RETURN 

5750 тар итә уру ив рһопес=хс=с=алпезлиелинлилтлдолтсни -"ылгаси 
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SEET EE SS 

270: begin: plot initialize əldəki. iii --...-............ 
9780 DISP "load plotter";@ INPUT WS ! dummy input 

ጋ ከኑ ከ 7 TEN 5320 G RETURN ! abort plot 

3800 PEN 1 @ FRAME € FENUP 

3810 RETURN 

5820 ! end: plot initialize 0፻0ር89ህ/8ወ።5።2555225222፡222፡፡፡6፡፡2::፡ 
5850 | begin: show plot Dprocedurgeszszzasszzzzzzzzzzzzzczzzzunzmncz 
5840 ! begin: scaling | 

кр иб 2,В5=-Б1 

39850 FOR 10=1 TO 15 | 

3870 WO=X (10) € Wl=Y(10) 

280 1F WO>B2 THEN 82=Ҹ0 

3870 IF WOXBI THEN BisWO 

5900 IF 4W1>B5 THEN 85=Ҹ| 

5910 IF W1«B4 THEN B4=54w1 

DOC NEXT 10 

aca range ¡bl do 


። 11 ከ. 01700 
3990 ! X axis scale 
39960 Bis5 ! w-interval 


3270 MOzS*B3 ! major divisions 
5:27 52 D TI 2/401-1P (51/7/0123 +45 
Seo if Sees THEN 5910 ELSE 5000 
ирс HEM 6080 ELSE 6160 
5010 ! too many divisions 

6020 K3=1 

5950 В5=В5*28(К5) @ KS=FNK (KS) 
5040 WO=5+B3 

X::7:755 кЕ1И82/И0)=-1Р(В1/Н0) 145 
Gos TREN 0) ELSE 5:00 
5070 ! end too many divisions 
5080 ! begin too few divisions 
8070 K3=7 

5100 B85:52/8(K3) 8 KO=FNI (KZ) 
71O5::55F5:: 60120 ELSE 6200 
5120 Н035482 

2777” El. 67 101-iP(21710172 
7:5775775:: - iHEİH 100 ELSE 6160 
6130 ! end too row divisions 

piso land if 
SID ANO 1 low «< biry 
7577577700 38 ! hi x bary 

50 50019955250 

7:1) method tailed=set detault x scaling 
25.0 5” 52-811710 

6220 88=10 
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6230 
6240 
5220 
6260 
6270 
6280 
5290 
5200 
5310 
632 
6330 
5340 
Gog 
6360 
6370 
2380 
6570 
2400 
5410 
5420 
420 
5440 
4450 
5450 
5470 
5480 
5490 
£300 
227 
6320 
SD 
8540 


= ^ 
x 


5360 
5570 
6530 
6370 
6510 
5610 
5529) 
6430 
5540 
5690 
6650 
6670 
5680 
5690 
6700 


A 


Cr 


DISP "default әс ии изи 

! end default scaling 

epnd x scaling 

| begin y scaling 

B6=.002 ! mV 

N2-7S*Bàó 

B9= (CEIL (23/42 I Obs o O 

LE 87:70 ተህ ሇፓሇ ክስ ም ምን፡፡፡፡ 

IF ከክ! የ ርክ፡ንሀ፡፡ ክህ: 

' begin too many divisions 

K3=3 | 

B6=Bö+B(k3) 8 42=5286 @ KS=FNK (KS) 
BEE EE 

İF 89740 THEM 5330 2158 01001 

! end too many divisions 

! begin too few divisions 

ሂ2= 

85-B5/B(KX3) G W2-5*Bà 8 K 
İF Bevo TREN 5420 FO Di 
EE S D EM 
Ir E EE ee 
! end too few divisions 
7:57 

84=1P (84/42) *Wž 
BS-EAB6£57 Tn YO) 
GOTO 6340 

) method fail-set default y-scaling 
Вӧ=(В2-В84)/10 

88-210 

A A или 

! end default y scaling 

! end y-scaling 

SCALE Ply bee oS 

!' end scaling 

1001010 610: 

EE 3110 

20810 7 

IF Koi! TREN kosy 8 0021033057 

FOR 10=1 TÜ Kə 

ከ DOCTR 1) 

NEXT 10 

UT RZ RBA 

KZ=K 3+ 

FOR Tosk TO TS 

PLOT 201 0 

NEXT TO 

БОТ 65 ии 

CSI EA O FAD 


20 


6710 
6720 
6730 
6740 
6750 
6760 
6770 
6780 
6790 
5800 


8580 


5810 
6820 
6830 
6840 
6850 
58606 
5870 
SEED 
5370 
6900 
6919 
550 
2950 
5740 
6950 
6960 
50 


505 a 1515 

N4=N4+1 8 H$z24$033,5371] 8 PRINT "plot ";H$;"R";N4 

SETGU 

Me DANOS OPA LORE 6 

LABEL USING "K" ; "amplitudelmV rms)" 

“ከከ ከ Li 4 © MOVE 72712 

IF S59=2 THEN W$="depthimm)" ELSE Vi-"Trequency(Hz)" 

LABEL USING UK! 3 KŻ 

MOVE 72 5 G SIZE ጋ 

ከ »ጹዱምኺፕቭጦኺጦኾጦሙሲፕቭዥ.ሲ0ንችፕ .: 55 00 000 8570 LUE 9$.,40$,N4 ELSE LABEL USING 
ЛОРА 

gc Tg Su ኤጌ» ጨጨ 

10:00 plot 

bean curve frt3plot 

Ел Bu 1 ÖÖ 

CUIDAD? STEP BS 

W2-*Wi-XO 

751177 Ur 7112770705 (12912 / (Ü7x09)) 

IF W2>B4 AND W3¢BS THEN PLOT W1,W3 

NEXT Wt 

FENUP 

111135 PEN 0 

36:02 €60f € Tit nlot 

7557127177 € cən. 010- 

12=0 G8 5920 


54:06 er 1161 
RE TURN 
мала" Bat рӱутосефигесс====алаатаапезтлаасллаллалттыагсас 


zı 


4730 
6970 
7000 
able 
7010 
7020 
7000 
7040 
OS 
7050 
7070 
7080 
7090 
7190 
7110 
2129 
m. 
7140 
7150 
7150 
500 
7 
7170 
0 
77 
7220 


ot 


ያለ 

BARD 
ZEW 
1280 
Ta TO 
7280 
7290 


7220 


| ======MEASUREMENT PROCEDURES==========================z=<= 
' begin: amplitude rəad setup procedures=s=s================= 
| input:56,F8,V6 out: timer #1 Set for resd enable, Pome 
Sl wt 

IF S6 THEN 7/0Z0TE CSES 3» 

begin) 2052105) 

! automatic drive adjust 

OUTPUT 707 USING 8800 ; V6,F6 İP 7201007066 E 

5120 8 ON TIMERE 1,T7*15*T7T9 GGGUB 7550. 70” 

! begin: loop to adjust VG until LIA output tee 
WAIT 2215 

ENTER 709 ş AS | read LIA magnitude-1gnore Fed m 

İF (57.5 THEN TOONE ESER TIO 

tLe Ue too həl 

IF AS+1.12 THEN Vó6=FNR(VE*. O) EGSE RV See he eee 

Bü” 

IF ASŚ.1 THEN 7150 ELSE 7290 ! hars s the 190) 6077707 
CTA OUT የሮሮ ә 
VAsMIN(ID,FNR(Yá¥15)) 
20197170 

| rasst HP5514 and retas 
8018010 70-45 (MGE 5 EI 
OFF TIMER# 1 8 51=0 

ON TIMERS 1,.,18+*(15- 54 5 6659 : 

PRINT Soret tt de" ede ive или 

| hold for read enable 

SIE NO US 111 

end: ADE block 
begin: no ADC block 
UTEU1 707 157110 ከ ከ ም ን 

120 8 ÜN TIMERE 1,T7+154715 G50SUE 72201 t ras p ри. 
епәтлә ALE ok 
<) 8 RETURN 
ənd: anplitude read setup procedure =።=፡፡፡222፡:222252መ፳2:-:=:2= 


- 
Va 


m 


—1 


A > A 


ያራ 


7310 
4520 
0 
7240 
2290 
7340 
7570 


1580, 


2290 
7500 
7410 
7420 
7430 
7440 
Arad 
7450) 
7470 
7480 
7496 
BANIE 
750 
noon 
gs 

7240 
5: 
7569 
157.0 
280 
Zo 0 
7690 


! begin:amplitude read proceduręzs====z======z============== 
Пипи WOSGUE V2, 3: 

! hold fcr read enable 

1500705 0 ELE 7550 

! read enable 

TRIBGER 723 ! measure xmtr voltage 

ENTER 709 ; AS ! LIA magnitude output 
V2=2A54*V5 ! volts receiver cutput 

EM VI mtr voltage 

A2=V2/V1 | 

! check for LIA overload 

İF “Aor to THEN 74S0MBESE 7450 

' alarm -overload 

7( 7U12717:15:2717127 PRINT FEAS: VT E BEEP 55,2000 
! end: overload test 

3155021757 ЕА 1011000 41 1000472 
507 107 NE 13450 TREN 7480 ELSE 72320 
graphics storage 

12=154| 

M pir 00-20 clst 11415)=Р5 
2775700002 ! mV 


! end graphics storage 

RETURN 

' end : amplitude read proceduressszzazsezzzuuzancussciomuscc 
| begin:raac enable procedurgss=zuizunzzzzzznunmumumizmwun mit 
' €alled by timerfi-resets read enablo flag sl 

ነ ከ1118. ! 


= Weg 
RETURN 
Ee El ersch Le proceduremssszaszzczczrsszszzzucuucmdsi 553 
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7610 ! begin: mova manipulator down 29 mmzzs2s2zzzzzzzzszz22zzzz2zz 
፣.ር ከ0 O 

7630 ENTER 720 ; NO ! counter en 

7540 15 0651070600 O SE 

7630 IF 2950 THEN OUTPUT 709 ¡"CLS1" Q Nos) ELE 007711 a ይዊ 
0" 8 N3s-1 8 29-020 058) 

7660 ON TIMERS 2,IP(29/R2-T2) B508UB 770007 52 RP TEE E 

71570 81,9220 Fredo on 

7580 RETURN 

7690 ! end: move дап1рајаҝог руссефигелх=елазлаалахданалпализаныын 
7700 ! begin: manipulator stop proceduresssssssssszzzoz2szm2:045 
Тилә 1024 eo WE ME 

7720 ! output revised 201 562205 

7730 IF N3>0 THEN QUIPUT 709 ¡"DPNT" ELSE ብህ ስግጥ ከ ን 
7720 DEF 1 h ди. 

ТӘЗӘ РЕКЛ бә op R 

//S0. 1000000705 ከክ 

7770 ENTS he cee 

7780 IF N2=N1 THEN. 7799 ELSE NI=N S m SEED 

E EEN 

7800 23:N3:(N1-N0Ö)/RÖ 

7810 20=20+28 

42:20 0140100017 ፡፤ 

2850 3270 MO o 

7340 * auto drive control 

¿230 IF Sa ከኮርዚሪረርጽቲ си ик 

7360 Pagin adju Orne 

{SO VENTER 705 S 

(SEU STF әз 2 HE 2570226 ጋ: 

PSG "Eg Балет 

7900 12 45>1,.12 THEN V6=—FNAUVG ©) 3) EECHER 

19 OR CUTE UT 7 E 

TI RIN o P 

1500 ти E 


YOLu TS 


79520 ON TIMER# 1,77+15+T5 GOSUB 7550 ! sst delay Lo “2210700 E 
70000 

777701 

7989 1 endimaninulatar ston 3rocedurasss=sz==s==========- ስከ 
¿970 * begin: manipulater relocate precadure======- 5.2 


EE 

01015 2050 THEN 00 
BOZU DISF "nove tarre INE EMI 
2060. 29=24-700l na 

GOZU IBOSUE 7610 1 move mm 

8050 54=: © RETURN 

8060 ! end: manipulator relocate procedure ии 


p 
t 
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8070 ! begin: title string justifiers====ss===s========re===9=->> 
8080 IO=LEN(AS) 
8090 IF IO MOD 32»0 THEN A$-A$&4" " @ GOTO 9080 
8100 RETURN 
81107 and title 41[1511፥125የ===2222222=፳2252222፳2፳222=22222:22:2:2=:22 7: 
ሸሽ. ADE Zn CE toggle A EE tur" 
EE EE 5170 
8140 ! operator call 
8120 $421 ! fast refresh menu 
27 0675021510 101 0 0 0 12 A 4=cal":2 INPUT. 19 
8170 ! process command 
571 0) 61766 sozuce Tds-DIg ! 19>0 turns ADC on 
eS MOD 5 
8200 IF 13=0 THEN 8420 
8210 OUTPUT 707 ¡"APOVO" 
8220 ON 18 8010 8230 3280 3330 8380 
2500757 
570 .00 020 ə ES LC 
s 51177 ss LS ING 882606 "EIA!" 
S SDSBUROUIC T USTND 0520: "9282" 
earn GORE 8420 
A 
кии RIMT 
75 0701177 5 UĞİNE 5820 + "B2A1" 
25777? 7... USTNE 5850 + "B2AZ" 
5-20 60710 8420 
2009101 
Ee eer TEEN 
777777 TUS INe 000 "ВСА2" ! 
7: 0115177711 USİNG 9820 : "ba" 
8270 GOTO 8420 
Bee? ' ofreame calibration 
би Se="F ical” 
1777170172: 7705 15146 0822 + Bal" 
GE E 714 USING 8820: "BTAL" 
Elle, "rf". Er 
77(( ::0 © RETURN 
9410 1 end: xmtr/adc toggle routinaesssezsszzzszuuzsnaunuuumuc- 
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8450 


8450 IMAGE /////K 


— IMAGER ити 
,50/'stepsi 


"ran 
8480 
8490 
8500 
Em 
8510 
8520 
8530 
8540 
3550 
8560 
2077 
Ди, 
3580 
70 


LO 
na 


5 


= 


E 
2 


tO (Er Du) m (1 — 


መሙ 
o™ 


— - 
c 


(Ə mM - 
О- O- ^ > i 
y U- “3 cn 

c 


сә 


a © 
> 


= 
z e 
— 


= 
መ” 

- 

መሥ 


-4 ir O~ 
coña 


CH 00 (à CO (3 0D CL 


~} -.4 
лә 


[3 <> 


pu 


(A vi OO 
= 

ud - 

Çə 


zd 


4D, 


get " 
IMAGE 
IMAGE 
IMAGE 
107 
IMAGE 
IMAGE 
IMAGE 
IMAGE 
IMAGE 
IMAGE 
IMAGE 
b. 
IMAGE 
IMAGE 
6 
IMAGE 
mments 
IMAGE 


IMAGE 
IMAGE 
/ 
IMAGE 
ALME 
[MAGE 
IMAGE 
IMAGE 
IMAGE 
IMAGE 
IMAGE 
IMAGE 
IMAGE 


,5D/"center 
50: 

ከ ዛዳ ሰ SIE 
4190420 
бәзи Ди 


“ҝә H 


>= € 


2D,"H2"/5X3D.5D,"nV"/5X8.3DE, 


/"resonance center +requene 01” 
ири 2 pt obscured E 

Ul os NO 

iə 7 2810 

"resonance quality factom BRE 
/'search mee Ee 
EE AE Ta 
/ 


OH 
= 


2. 


eene 
D, EE 


DAD DA 
[UU REI Ы 
TU RUE mE RUE 
/"üsreturn to main 

-tonstants" 
"lockinsansl67 7 ከ0 
OSOBA 0 
"depth (na) :" 
5717100106”) 


ግ 


Tee 


gə” 


12[/5 ኬክ 


000 4D, "mV" 7 


25 m 


"volume 
bo c. 


ግጦ” 


^ 


11 m 


LEE 
AD 


dee = 


BAM 
xd: 

“m MET 
2 6 1 Ge 
1507 
ጣጫ 1 
ә, 
K PUT 
DE 


ገን aa, 


ONUS 
EFEKTE", 
8" SIE 

ptz"/"spe=",4D,2D,"+/-", 
corra", s pm: 


8 4DE) 


መዉ 
0.20.20. 


SDr በብር c mi sm cc 


" 


"niaZnasx "dime cd rm 


, 4D. D, 


it "1 Lé 


gö DN 
mus Eege 
20 spd, 23 t 
lera EAD ` 
i NB AH 7 s: 


preamp 
p AES E 
7 
була ,K/"dens 
Tp 45.20," dB re 
5: ES 2 A e 
АР" 28 5) 11 
"AFA ORA 


y 
V ар 


иә 
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treg measurement"/"estimate:” 


ирили 
2 ) 


Y 
* 


a a 9 4^ 4 ao yh መወሠ ፦ሠወ መሕ መህ vm 
. 


! Images for formatted outputszssszzzzzzzzzzzzzzcnn 


| 9D 


"nV /H2^ 2 M 


re ; 
J d: ብ. . 


qi dn 
. SE 


“ıı 


O 


L 
^ 


-” wom pe w 


- 


` = U 


ә u UP a q 4 / 


8320 IMAGE %,4A 
3330 | end: 1 0326555522ሪ2222222222222222222222222=:22፡2 225፳2=:22=:: 
8840 ! begin: end program routinezsssszszzzzzz225zz25zz5zi:z5zzzzrze 
230 0UTRUT 709 4 "RS" 

: ከ ENER? 

3870 END 
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